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ABSTRACT 

We present 7 spectroscopically confirmed Type II cluster supernovae (SNell) discovered in the 
Multi-Epoch Nearby Cluster Survey, a supernova survey targeting 57 low redshift 0.05 < z < 0.15 
galaxy clusters with the Canada-France-Hawaii Telescope. We find the rate of Type II supernovae 
within i?2oo of z ~ 0.1 galaxy clusters to be 0.026t°;oig(stat)to;°°i(sys) SNuM. Surprisingly, one 
SN II is in a red sequence host galaxy that shows no clear evidence of recent star formation. This is 
unambiguous evidence in support of ongoing, low-level star formation in at least some cluster elliptical 
galaxies, and illustrates that galaxies that appear to be quiescent cannot be assumed to host only Type 
la SNe. Based on this single SNII we make the first measurement of the SNII rate in red sequence 
galaxies, and find it to be 0.007+o'oo7( s t a t)lo'ooi( s y s ) SNuM. We also make the first derivation of 
cluster specific star formation rates (sSFR) from cluster SNII rates. We find that for all galaxy types, 
the sSFR is 5.ll3 5 1 8 (stat)±0.9(sys) M Q yr" 1 (lO 12 M )-\ and for red sequence galaxies only, it is 

2.0;o;g(stat)±0.4(sys) M yr" 1 (10 12 M Q ) _1 . These values agree with SFRs measured from infrared 
and ultraviolet photometry, and Ha emission from optical spectroscopy. Additionally, we use the SFR 
derived from our SNII rate to show that although a small fraction of cluster Type la SNe may originate 
in the young stellar population and experience a short delay time, these results do not preclude the 
use of cluster SN la rates to derive the late-time delay time distribution for SNe la. 
Subject headings: supernovae: general — galaxies: clusters 



1. INTRODUCTION 

The rates and properties of Type II supernovae 
(SNe II) indicate they are explosions induced by the col- 
lapse of iron cores in stars of initial masses 8 < M < 
20 M (|Smartt et alj 12009). Such massive progenitors 
have been directl y confirmed in pre-explosion images 
for several SNe II (lLi~eTal.ll2007al [Gal-Yam et all 120071: 
Elia s-Rosa et al. 11201 lUMaund et al.l 
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20111 ). Stars of initial mass > 8 M Q explode as SNe with 
a delay time, the time between star formation and explo- 
sion, of < 30 Myr (e.g. Henyey et al. 1959). Due to this 
relatively short delay time the SNII rate, SNRjj, is a 
direct indication of the current star formation rate, SFR 
(Botticella et al. 2012). Since SNe II are bright, using 
SNRjj to identify very low levels of star formation (SF) 
can be advantageous to using ultraviolet (UV) photome- 
try or optical spectroscopy. For example, small amounts 
of UV light may be undetectable in a luminous ellipti- 
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cal, and a little Ha emission in an elliptical can be over- 
whelmed by an active galactic nucleus. Also, at low red- 
shift UV photometry can be an ambiguous SF tracer be- 
cause of contributions from blue horizontal branch stars 
(the UV upturn) . The SNRjj is thus an especially valu- 
able SFR proxy in rich galaxy clusters, which are mainly 
composed of luminous elliptical galaxies. 

Evidence for low levels of SF has recently been 
detected in low redshi f t clu ster galaxies from opti- 
cal spectra (iFritz et al.l I2011P) and infrared photome- 
try ([Chung et al. 1 1201 If ) , and in low red shift field red 
seque nce galaxies from UV photometry (|Kavirai et al.l 
I2010T) . Despite this, SNe II have rarely been observed in 
early-type galaxies wher e the bulk of the stellar mass is in 
old stellar populations (Hakobv an et al.l l2008). As a re- 
sult, the SNRu in field ellipticals and cluster galaxies is 
not well constrained. Table [T] presents the current litera- 
ture values for SNRu. The Lick Observatory Supernova 
Search (LOSS; Leaman et al. 2011) placed an upper limit 
on SNRu in field ellipticals, and provi ded the first mea - 
surement of SNR H in field SO galaxies (|Li et aLll20TTH ). 
Five SN surveys were combined by Mannucci et al. (2008; 
hereafter M08), who measured the SNRu in galaxy clus- 
ters to be half the rate in field galaxies, but well above 
the upper limits of field ellipticals. However, the surveys 
complied by M08 were biased towards the most massive 
cluster members, and their SNe were not all spectroscop- 
ically classified. In this paper we present the SNRjj 
measurement from our large, complete, well character- 
ized cluster SN survey at low redshifts, and make the 
first comparison of SNII-derived cluster SFR to the re- 
cently detected low levels of SF from optical, IR, and UV 
data. 
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In contrast to SNell, SNela occur in both young 
and old stellar population s ()Mannucci et al.l 120051 : 
iScannapieco fc Bildstenll2005D . On going star formation 
in galaxy clusters is a concern when using the cluster 
SNRj a to constrain the slope of the SNIa delay time 
distribution (DTD). The two leading scenarios for the 
SNIa progenitor are a carbon-oxygen white dwarf ac- 
creting material from a main sequence or red giant star 
(single degenerate), or accreting from or merging with 
another white dwarf (double degenerate). Each occur 
over different timescales and predict distinctive DTDs at 
late times; i.e. the double degenerate scenario predicts 
more delayed SNe la. The colors of cluster red sequence 
galaxies indicate that their star formation was truncated 
at high redshift, and they have evolved passively since 
(e.g. Stanford et al. 1998; Eisenhardt et al. 2008). 
Based on this, it is assumed that cluster SNela have 
all experienced long delay times, and that the cluster 
SNRj a (z) can constrain the late-time DTD (i.e. Maoz 
et al. 2010). However, the presence of ongoing star for- 
mation suggests the cluster population may be contam- 
inated by short-delay SNe la. In this paper we use our 
SN Il-derived cluster SFR to evaluate this possibility. 

The Multi-Epoch Nearby Cluster Survey (MENeaCS) 
surveyed 57 low redshift 0.05 < z < 0.15 galaxy clusters 
for two years with the Canada-France-Hawaii Telescope 
(CFHT). Within the virial radius, -R200, of our cluster 
sample, we spectroscopically classified 7 cluster SNe II 
and 23 cluster Type la supernovae. This paper is part of 
a series based on MENeaCS. In Sand et al. (2011; here- 
after Sll) we use the relative number of hosted and host- 
less SNe la to determine the mass fraction of intracluster 
stars. In Sand et al. (2012; hereafter S12), we mea- 
sure the cluster SNRj a from MENeaCS, and combine 
it with published SNIa rates between 0.02 < z < 1.12 
to constrain the slope of the SNIa delay-time distribu- 
tion. Two additional MENeaCS papers are nearing pub- 
lication: one showing evolution in the cluster dwarf-to- 
giant galaxy ratio over redshift (Bildfell et al. 2012), and 
one investigating the demographics of tidally disturbed 
galaxies in clusters (Adams et al. in preparation). 

In this work we present the 7 cluster SNe II, including 
the unprecedented occurrence of a SN II in a red sequence 
galaxy. In § [5] we describe MENeaCS, and present the 
photometric and spectroscopic properties of our 7 cluster 
SNe II and their hosts galaxies, including a comparison 
of SNII and SNIa hosts. In §[3] we constrain the level of 
ongoing star formation in our one red sequence SN II host 
galaxy using published multi- wavelength data. In §@]we 
describe our SNII rate calculation and its uncertainties, 
and compare our results to published rates. In §[5] we de- 
rive the star formation rate in clusters, and compare it to 
measurements from IR, UV, and spectroscopic observa- 
tions of cluster galaxies. We also discuss the implications 
for the SNIa DTD. We provide a summary of this paper 
in § |51 and in all cases we use a standard flat cosmology 
of (0 M , fi A ) = (0.3, 0.7) and h = 0.7. 

2. THE MENEACS SN II SAMPLE 

MENeaCS monitored 57 rich galaxy clusters with 
monthly cadence for two years in the g'- and r'-band fil- 
ters with MegaCam (|Boulade et alj |2003) at the CFHT. 
The survey design, cluster sample, observing strategy, 
real-time reductions, transient detection pipeline, photo- 



metric calibration, host identification, and spectroscopic 
SN classification techniques we used to discover super- 
novae are presented in S12. The data were reduced and 
searched for transients in real time, with regular spec- 
troscopic runs scheduled each month of the survey for 
supernova classification with the Blue Channel Spectro- 
graph (BCS; Schmidt et al. 1989) or Hectospec (Fabri- 
cant et al. 2005) at the MMT Observatory. Photometric 
calibrations to SDSS filters g and r in the AB magnitude 
system were performed using standard stars in the fields. 

We spectroscopically followed up all SN candidates 
brighter than m g =22.5 magnitudes, and with colors 
g — r > 0.8. As described in S12, we used the pub- 
licly available Supernova Identification (SNID) routine 
of Blondin & Tonry (2007) to spectroscopically classify 
our SNe. Cluster membership was assigned for SNe with 
\vsn — v c i uster \ < 3000 km s" 1 . In total we confirmed 
23 SNela (4 of which were hostless), 7 SNe II associ- 
ated with our clusters, and 37 background SNe. Due to 
the one square degree field of view of MegaCam, we are 
complete to ~ -R200 and can calculate SN rates within 
this radius. For comparison with past surveys, we also 
present rates within R = 1 Mpc. 

Table [2] presents the 7 SNe II discovered in MENeaCS 
clusters including their internal identification name; the 
UT date of spectroscopy; the telescope and instrument 
for follow-up; cluster redshift; galaxy redshift where 
available; details of the SNID best fit including red- 
shift (and uncertainty) , SN template type and name, and 
median phase (and the standard deviation); and finally 
the spectroscopic exposure time in seconds. Figure Q] 
presents the classification spectra for each; all 7 are best 
fit with Type II plateau spectral templates. 

As discussed by Li et al. (2011a), reliable distinction 
between Type II subtypes requires multiple spectra and 
well sampled light curves. Type Iln (narrow spectral 
lines) are spectrally distinctive but can evolve to resem- 
ble a regular Type II. The Type ILL are spectroscopi- 
cally similar to Type IIP, but have a linearly declining 
light curve and no plateau phase. Type lib (broad spec- 
tral lines) can resemble normal Type II at early times, 
but their light curves are distinctly double-peaked. With 
the single epoch of spectroscopy and monthly photomet- 
ric cadence of MENeaCS, we cannot confidently identify 
subtypes for our SNe II. 

2.1. MENeaCS SNII Host Galaxies 

Photometry for MENeaCS SNII host galaxies is mea- 
sured from SN-free deep stack images, as described in 
S12. Table |3] presents the host galaxy details for each 
cluster SN II, including the coordinates, r-band magni- 
tude, g — r color, the A(g — r)ns color offset from the 
cluster's red sequence, and the clustercentric radius in 
units of kiloparsecs and i?200- Clustercentric radius is 
the distance from the brightest cluster galaxy (BCG); 
for our 7 SNe II host clusters the BCGs are within 30" of 
the X-ray centers (from Chandra and ROSAT), except 
for Abell 2443 which has a ~ 1.5' offset. However, even 
a potential 1.5' shift could not cause any SNe II to cross 
the R200 boundary, and be included or excluded from the 
sample. 

Given that SNII progenitors are young, but SNela oc- 
cur in both young and old stellar populations, we ex- 
pect the distributions of g — r colors to be different for 
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SNII and SNIa host galaxies. In Figure [5] we plot the 
A(g — r)us color offset from the cluster red sequence as a 
function of r-band magnitude for SNIa and SNII hosts. 
The dashed lines represent the median scatter in the red 
sequence for the MENeaCS sample. All host galaxies 
with A(g — r)ns error bars overlapping this zone are 
considered to lie "on" the red sequence, and all others 
lie "off' the red sequence. Unlike the SNela host pop- 
ulation, SNII hosts lie off their clusters' red sequences 
- except for the host SNII Abell399_ll_19_0, discussed 
in § El In Figure [31 we show that the distributions of 
host g — r colors are significantly different for SNe la and 
SNe II. The KS-test probability that the two samples are 
drawn from the same underlying distribution is low, but 
not negligible, at ~ 8%. 

Chung et al. (2011) show that the fraction of star form- 
ing galaxies increases with projected clustercentric radius 
within i?200- We therefore expect the radial distribution 
of SNII hosts to differ from SNIa hosts, which should 
in turn follow the cluster luminosity profile. In Figure [2] 
we plot the A(g — r)ns color offset from the cluster red 
sequence as a function of projected clustercentric radius, 
for both SNIa and SNII hosts. We note that no off-RS 
MENeaCS hosts are observed within 0.3 i?200> consistent 
with the view that the star formatio n fraction decrease s 
in the central regions of clusters (jChung et al. |[20TT1 ). 
Figure |3j shows the distributions of projected cluster- 
centric radii for SNIa and SNII hosts. The normalized, 
cumulative fraction of cluster g-band luminosity in red 
sequence galaxies is also plotted. Within -R200 it appears 
the number of SNe is roughly proportional to luminosity. 

3. THE RED SEQUENCE SNII HOST 

Unexpectedly, one of our seven MENeaCS SNe II oc- 
curred in a red sequence galaxy. As we discuss below, 
based on its spectrum and photometry we arc confident 
that it was a SNII and not a SNIa. The color, magni- 
tude, red sequence offset, and clustercentric distance for 
the host galaxy of Abell399.11_19.0 are listed in Tabled 
and are shown to be consistent with the red sequence in 
Figure [2] In Figure [4] we show an image of this galaxy 
with isophotal contours to highlight this galaxy's ellipti- 
cal morphology. 

When core collapse SN (CCSNe, Types II and Ibc) 
are discovered in elliptical galaxies, further inspection 
almost always reveals the presence of star formation. 
Hakobyan et al. (2008) appraised 22 elliptical galaxies 
hosting CC SNe; of them, 19 were misclassified as ellipti- 
cal, and three showed evidence of mergers or interactions 
and were thus likely to harbor recent star formation. Suh 
et al. (2011) investigated the near- ultraviolet and radio 
properties of nine early- type CC SN hosts, finding clear 
evidence of recent star formation in all. These results 
are consistent with the relative youth of SN II progenitor 
stars. In this section we look for evidence of recent star 
formation in the UV, IR, spectral, and radio properties 
of the host of Abell399_ll_19_0. 

3.1. SN Classification 

The classification spectrum taken for Abell399_ll_19_0 
with Hcctospec at the MMT Observatory is shown in Fig- 
ure [3 along with the best fitting SN template spectrum 
from SNID (|Blondin fc Tonrvl [20071 ) . The P-Cy gni pro- 
file, distinctive for SNII, appears at ~ 7000 A. Since 



the fit with SNID may not be overwhelmingly convinc- 
ing for all our readers, we do an additional analysis with 
Superfit ([Howell et al.l [20051 ). This routine achieves a 
better looking fit because it removes the host galaxy 
spectrum. First we run Superfit with loose redshift con- 
straints, 0.05 < z < 0.25, to independently confirm the 
SN redshift. We know the host's spectroscopic redshift is 
z = 0.072, and the SN-host association is unambiguous 
(Figure [4]). The top 3 best fits are SNe II at z = 0.07. 
We then run Superfit two additional times with the red- 
shift constrained to z = 0.07: first allowing only SNII 
templates, then SNIa only. The best fits, in Figure [SJ 
show that a SN II is the better match. 

The MENeaCS cadence of one epoch per month did 
not generate well sampled light curves, but given the 
relative importance of Abell399_ll_19_0 we discuss its 
photometry briefly Following a non-detection epoch in 
October 2009, this transient was detected in three con- 
secutive months, after which the MENeaCS observations 
ended. This transient had m g ~ m r ~ 20 magnitudes 
in the detection epoch (corresponding to M ~ —17.5 
magnitudes) . In ~ 55 days it declined by Am g ~ 2 mag- 
nitudes and Am r ~ 1 magnitude. When compared to 
the light curve templates of Nugent et al. (2002) and 
SN luminosity functions of Li et al. (2011a), this color, 
magnitude, and slow decline are all most consistent with 
a SNII Plateau. Regular SNela are too bright, and faint 
SNela decline too quickly (Phillips 1993; Perlmutter et 
al. 1997). Finally, the preceding non-detection epoch 
prohibits this from being the late-time shallow-decline 
epochs of a SN la. 

3.2. Ultraviolet Photometry 

If this host galaxy experienced a small amount of recent 
star formation, it might be evident in its near-ultraviolet 
photometry. Schawinski (2009) generate model spec- 
tra from near ultraviolet (NUV) to optical wavelengths 
by parametrizing the star formation history of early- 
type galaxies as a large population of old stars plus a 
small amount of young stars. They show how a galaxy's 
NUV — r color indicates the time elapsed since the most 
recent burst of star formation, given the fraction of stellar 
mass synthesized in the burst. For example, in an early- 
type galaxy where ~ 1% of the stellar mass is < 50 Myr 
old, NUV — r < 1. Similarly, Kaviraj (2010) use star for- 
mation history models to derive a relationship between 
the UV and optical photometry of bulge-dominated red 
galaxies, and the age and mass fraction of their most 
recent epoch of SF. 

We look for the UV counterpart of this host in the 
GALEX data release £3- There is no coincident ob- 
ject in the catalog, and under visual inspection the tiles 
show no hint of a source. From Bianchi et al. (2011), 
we know that UV sources are detected at 5a down to 
NUV ~ 20.8 magnitudes in the GALEX All-sky Imag- 
ing Survey, and that the Medium-depth Imaging Survey 
does not cover this region of sky. With this limit, wc 
restrict the NUV — r color to > 4.5. Based on Schaw- 
inski (2009), this constrains the fraction of stars < 50 
Myr old to < 1% of the total stellar mass of this SNII 
host galaxy. Comparably, the work presented in Kaviraj 

10 http://galex.stsci.edu/GR6/ 
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(2010) constrains the age and mass fraction of the most 
recent burst to > 250 Myr and < 0.5%, respectively. 

3.3. Infrared Photometry 

Photometry at optical and near-infrared (NIR) wave- 
lengths can be combined to reveal recent star formation 
in an evolved galaxy. Stellar population synthesis mod- 
els have been used by Li et al. (2007) to determine that 
B — V and B — K photometric colors best disentangle 
the degeneracy between galaxy age and metallicity. They 
show that a galaxy's location on a B — V vs. B — K 
color-color plot can be used to estimate the fraction of 
mass in young stars. To obtain these colors for this 
host we begin with our photometry, m g = 17.1 ± 0.02 
magnitudes and color g — r = 0.85 ± 0.04, and add 
2MASS m K = 13.62 ± 0.18 magnitudes (from NASA 
Extragalactic Database). We then K-correct to z = 
(jChilingarian et all I2010T ) . and a pply filter transfor ma- 
tions derived for stars at z = (| Jester et al.l 120051) . to 
obtain B - V = 0.9 ± 0.04 and B - K = 3.6 ± 0.2. These 
colors su ggest < 0.5% of the stellar mass is younger than 
0.5 Gyr (|Li et "alJl2007bh . 

Galaxies which appear quiescent from optical and NIR 
photometry can harbor dust-obscured star formation, 
and be luminous at far-infrared (FIR) wavelengths. In 
this scenario, dust absorbs UV light and re-emits it in 
the FIR. Chary & Elbaz (2001) present conversions from 
the Spitzer-MIPS 24 /xm wavebands to infrared lumi- 
nosity, LiR, & n d star formation rates, SFRir. For 
example, Graham et al. (2010) applied these conver- 
sions to Spitzer-MIPS fluxes for 20 optically elliptical 
SNIa host galaxies from the Supernova Legacy Survey, 
and found 2 were actually Luminous Infrared Galaxies 
(LIRGs, Ljp> > 10 11 L©) with specific star formation 
rates ~ 500 M Q yr" 1 (IO^Mq)^ 1 . 

This red sequence host was observed by the Wide-field 
Infrared Survey Explorer (WISE; Wright et al. 2010) in 
its mission to create an all-sky infrared map, and the 
WISE-W4 filter at 22 nra compares well with Spitzer 
MIPS 24 /in0. In W4, this galaxy has an apparent 
magnitude of 15.2 in the AB system (8.6 in the Vega 
system), but a signal-to- noise ratio of just 1.4 and a a = 
null, indicating this magnitude is a 95% confidence upper 
limit. We convert this magnitude limit to a flux and find 
Lin < 10 9 L Q . This is not a LIRG masquerading as a 
quiescent elliptical. 

3.4. Optical Spectroscopy 

We obtained an optical spectrum of this galaxy with 
the Blue Channel Spectrograph at the MMT Observa- 
tory as part of our program to gather spectra for all our 
SNela cluster hosts. Full spectral analyses will be per- 
formed for all of our SN cluster host galaxies in future 
work (Graham et al. 2012, in preparation). 

The partial spectrum presented in Figure |4] reveals 
hydrogen emission, indicative of star formation - but 
also shows the nitrogen, oxygen, and sulfur signa- 
tures of a low-ionization nuclear emission-line region 
(LINER). Kewley et al. (2006) show that star for- 
mation is not the dominant source of emission when 
\og{[NII]/Ha) > 0.0 and \og{[S 1 1] / H a) > 0.0. A sim- 
ple analysis of this galaxy's line intensities finds that [Nil] 

11 Supplement at http://wise2.ipac.caltech.edu/docs/release/prclim 



and [SII] are stronger than Ha (log([NII]/Ha) ~ 1.1 
and log([SII]/Ha) ~ 0.6). We estimate the maximum 
amount of Ha absorption by fitting template spectra of 
elliptical galaxies (Kinney et al. (1996); Fioc & Rocca- 
Volmerange 1997). After accounting for the template fit 
with the largest absorption, we are confident that the 
line intensity of Ha is < 4. This indicates the min- 
imum line ratio values are log ([NI I] /Ha) > 0.4 and 
log([iS I 1] I H a) > 0.0, which is still consistent with a 
LINER. While these spectral emission lines mean we can- 
not attribute the Ha to SF, we also cannot rule it out. 
The slit did cross the galaxy core, and future observa- 
tions with a slit orientation avoiding the core may reveal 
SF at this SN's location. 

3.5. Radio Power 

Radio-loud emission (Li.4ghz > 10 29 erg s _1 Hz -1 ) 
from elliptical galaxies could indirectly represent ongoing 
star formation (e.g. Delia Valle et al. 2005). We checked 
published radio source catalogs for 1.4 GHz emission at 
this galaxy's position. Its coordinates are not covered 
by the VLA FIRST SurvejQand it wa s not detected in 
the NRAO VLA Sky Surve£3 (NVSS). The complete- 
ness limit of the NVSS at 1.4 GHz is 2.5 mJy. At the 
rcdshift of this host, z — 0.072, the radio-loud popu- 
lation is incomplete and we cannot constrain the radio 
properties of this host. 

In summary we find no evidence of star formation in 
this galaxy, aside from the presence of the SNII. This 
indicates that either very low levels of star formation 
and trace amounts of young stellar populations can exist 
in red sequence galaxies, or there is a rare other channel 
to SNe II with a longer delay time. 

4. THE CLUSTER SN II RATE 

Our calculation of the SN II rate in clusters follows the 
method we used for SNIa rates in S12, which is very 
similar to that used for high and low redshift SN rates 
by Sharon et al. (2007) and Barbary et al. (2012). The 
rate of Type II supernova, SNRjj, is calculated by: 



SNR n = 



Nn x C mc /C ; 



spec 



(1) 



where Njj is the observed number of SNe II. The spectro- 
scopic completeness, C spec = 0.91, accounts for the ~ 9% 
of the time when MENeaCS was detecting SNe, but we 
did not have spectroscopic follow-up due to weather and 
telescope scheduling. This value is independent of SN 
type (see Sll and S12). The inclination correction, Ci nc , 
accounts for SNe II that are undetectable due to extreme 
dust obscuration in highly inclined and edge-on spiral 
galaxies (e.g. Cappellaro et al. 1993b; Cappellaro et 
al. 1999); our inclination correction, d nc — 1.62, is de- 
rived in Appendix A. Over all N ep survey epochs of every 
cluster we sum the control time for that epoch, Atj , mul- 
tiplied by the mass or luminosity surveyed in that epoch, 
Mj (see S12 for a description of how these are calculated 
from our deep image stacks). 

The control time is the effective amount of time sur- 
veyed by the j th epoch, expressed by: 

12 http://sundog.stsci.edu/ 

13 http://www.cv.nrao.edu/nvss/ 
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At = / r)(m(t))dt, (2) 
Jti 

where rj(m), the MENcaCS detection efficiency as a func- 
tion of apparent magnitude, is determined from simu- 
lated transients and is presented in S12. Although the 
population of simulated transients used for our recov- 
ery statistics have magnitude distributions that mimic 
a sample of SNela, the resulting detection efficiency is 
appropriate for use with our SNell. For the SNII light 
curve, m(t), we start with the absolute V-band SNII 
template light curves from Nugent et al. (2002). ME- 
NcaCS detection efficiencies were calculated in the g- 
band. We convert from My to m g using the cluster's 
redshift, the SN II K-correction (based on spectral tem- 
plates from Nugent et al. 2002), and the photometric 
calibrations for the j th epoch. The integration bound- 
aries t\ and ti are defined by the time during which the 
SN II template light curve meets our color limit for spec- 
troscopic follow-up, g — r < 0.8 (e.g. 39 days for z = 0.15, 
and 66 days for z = 0.05). We account for potential re- 
discoveries of the same transient in multiple epochs by 
subtracting from r](m) the probability that it was de- 
tected previously, in the same fashion as SI 2 and Sharon 
et al. (2007). 

We use a Monte Carlo method in which the rate is cal- 
culated many times. For each realization a peak absolute 
magnitude is randomly chosen from the luminosity func- 
tion discussed in § 14.11 At every instance of Mj, r](m), 
and Nii, we randomly draw their value from an appropri- 
ate distribution based on their uncertainty (e.g. Poisson 
error for Nn). We run this Monte Carlo calculation for 
500 realizations, which generates a distribution of rates. 
The final value for the rate is the median of this distribu- 
tion, and the statistical uncertainties correspond to the 
16th and 84th percentiles (the 68% confidence interval). 
The final results are presented in § 14.31 

4.1. The SN II Luminosity Function 

For each realization of the Monte Carlo we randomly 
draw the SNII subtype (P, L, b, or n) and peak absolute 
magnitude from the volume-limited luminosity functions 
(LFs) published by the LOSS (Li et al. 2011a; hereafter 
Lilla). We use their LFs for type SO-Sbc hosts because 
the majority of the stellar mass in cluster environments 
is in early-type galaxies. The fractions of each SNII sub- 
type in SObc host galaxies are: P, 72%; L, 12%; b, 9%; 
and n, 7%. These fractions do not change by more than 
2% when all host types are considered, and the SNII 
LF for all host types is very similar to that for SObc 
only. We also use the appropriate light curve template 
for each subtype from Nugent et al. (2002); the SNIIL 
light curve when subtype L is chosen, and the SNIIP 
light curve when subtypes P, b, or n are chosen. 

To compare to the SNRn from M08, we repeat the 
Monte Carlo process using the same SNII LF as them: 
that of Cappellaro et al. (1993a; hereafter C93a). They 
present the SNII LF as Gaussian functions for sub- 
types P and L separately: peak Mb,iip = —16.38, 
anp = 1.49, peak Mb.iil = —16.82, and am = 1.1 
magnitudes. In Figure [5] we compare the SNII LFs 
from Lilla and C93a. To plot a single LF from C93a, 
we combine the IIP and IIL into one Gaussian of peak 



Mb = —16.51 magnitudes, a — 1.85 magnitudes, assum- 
ing 30% SNIIL and 70% SNIIP. As is evident in Figure 
[H Lilla detects a population of faint SN II-P (M s ~ -14 
magnitudes), which results in a non-Gaussian LF. How- 
ever, since MENeaCS is not sensitive to transients fainter 
than Mb ~ —14.6 magnitudes, these two LFs produce 
effectively similar results. This is discussed further in 

Here we make two important notes about how we in- 
corporate the SN II LFs and light curves. First, the Lilla 
distribution of absolute peak magnitudes from their un- 
filtered survey is very closely matched to R-band, and 
can be considered as Mr magnitudes with no correction. 
The C93a LF is for Mb, which peaks several days earlier, 
and the Nugent et al. (2002) light curve templates are 
for My . Fortunately, we do not need to convert between 
filters because the intrinsic B — V and V — R colors of 
SNe II at the time of B- and V-band maximum light is 
~ (e.g. Poznanski et al. 2002; d'Andrea et al. 2010). 
Although SNe II will be redder a few days later, at the 
time of maximum light in the i?-band this is due to a de- 
cline in B and V; the i?-band magnitude increases only 
slightly between the times of B- and i?-band maximum 
light. Therefore, we directly apply the i?-band LF of 
Lilla and the S-band LF of C93a to the V-band light 
curve templates of Nugent et al. (2002). 

Second, neither Lilla nor C93a correct their SNII 
LFs for host extinction, and by choosing randomly from 
these LFs we automatically include host extinction in 
our Monte Carlo rate calculation. By using the Lilla 
LF for SO-Sbc type hosts, the host extinctions are as 
similar to that expected for cluster galaxies as possible 
because most cluster galaxies are of similar early types. 
While the observed colors of SNe II-P do have a spread 
due to host reddening (e.g. Hamuy 2003; Krisciunas et 
al. 2009; Olivares E. et al. 2010), this affect is expected 
to be small for most of our surveyed mass in galaxy clus- 
ters. For example, Hamuy (2003) find that SNe II asso- 
ciated with groups/clusters show little to no reddening. 
The affect of host reddening and extinction on SNRu 
is accounted for by the combination of our chosen LFs 
because they are uncorrected for host dust, and by our 
inclination correction which is discussed below. 

4.2. Interlopers 

As mentioned in § [2J cluster membership was assigned 
for SNe with \v$n — ^clusterl < 3000 km s -1 , which ac- 
tually includes ~ 50 Mpc in front of and behind each 
cluster. Any SNe II exploding in Hubble flow galaxies 
within this cylindrical volume may be erroneously associ- 
ated with our galaxy clusters. The number of interlopers 
we expect to have observed, N exp , after accounting for 
our detection efficiencies is: 

r r 3=N ' P 

Nexp = vo1 s P ec £ AtjVj _ (3) 

J = l 

The volumetric SN II rate at z ~ 0.1 is R vo i ~7±3x 
10~ 5 SN yr- 1 Mpc" 3 (jBazin et al.l2009fl . The MENeaCS 
spectroscopic completeness term C spec = 0.91. The in- 
clination correction, C, nc , is discussed in Appendix A. 
We assume an interloper-hosting field galaxy would not 
be elliptical, which slightly raises the inclination correc- 
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tion factor used here to Ci nc = 1.72. The final term is 
the control time, Ai, multiplied by Vj, the comoving vol- 
ume element within ±3000 km s^ 1 and the chosen cluster 
radius (1 Mpc or i?20o)- This term is summed over all 
observed epochs, N ep . Our control times are shorter than 
the 1 year of MENeaCS survey time because SNe II are 
intrinsically fainter than the SNe la which MENeaCS was 
designed to find. 

For a cluster radius of 1 Mpc, the result is an ex- 
pected number of interloping supernovae N exp — 0.2j^'^. 
The Poisson probability that we observed 0, 1, or 2 SNe 
within 1 Mpc is ~ 0.81, 0.17, and 0.02 respectively. Sim- 
ilarly for a cluster radius of -R200, N exp = 0.7l ; 3 , and 
the probabilities of observing 0, 1, 2, or 3 SNe are ~ 0.50, 
0.35, 0.12, and 0.03. For every realization of our Monte 
Carlo we randomly draw a value of N exp = 0, 1, 2, or 3 
SNe, weighted by its respective probability, and subtract 
it from the number observed in order to produce SNII 
rates statistically corrected for interloping SNe. Since 
our one red sequence SN II host has a redshift consistent 
with the cluster and lies on the photometric red sequence 
- and considering that the number density of potential 
red sequence interlopers is relatively small - we consider 
it very unlikely to be an interloper, and we do not apply 
the interloper correction to the red sequence SNRji. 

4.3. MENeaCS SNII Rates and Uncertainties 

The MENeaCS SN II cluster rates are presented in Ta- 
bleful the conventional units of SNuB and SNuM, where 
SNuB = SNe(100 yr 1Q W L B ^)- 1 and SNuM = SNe(100 
yr 10 10 l Af Q )~ 1 . We present the rate in three cluster 
galaxy subsets: "All", the total stellar mass including 
the intracluster stars; "RS" , red sequence galaxies only; 
and "Off RS", galaxies lying off their cluster's red se- 
quence. We show our results with the SNII luminosity 
functions from both Lilla and C93a, and within cluster 
radii of 1 Mpc and i?200- 

The statistical uncertainties include the Poisson error 
on the number of SNe II observed, the interloper con- 
tamination, the uncertainty on cluster luminosity, and 
the uncertainty in our detection efficiencies. Statistical 
errors are dominated by the Poisson uncertainties and 
interloper contamination (with relative contributions of 
~2/3 and ~l/3, respectively), with very small uncer- 
tainties (1~5%) from our detection efficiency and cluster 
mass/luminosity uncertainty. The systematic uncertain- 
ties include a small contribution from the inclination cor- 
rection factor (5%; Appendix A), an offset of ~ 10% in 
cluster luminosity which was derived from a comparison 
of MENeaCS photometry to SDSS, and the uncertainty 
in i?2oo for rates within i?2oo (the latter two are discussed 
in detail in S12). The relative contributions to the sys- 
tematic uncertainty from these two components is ~ 30% 
and 70% respectively. 

4.4. Caveats 

Here we make several important notes regarding the 
MENeaCS SNII rate in clusters. First, ZwC10628_7_08_0 
was both m g > 22.5 magnitudes and g — r > 0.8 at all de- 
tections, and therefore did not officially meet our follow- 
up requirements. Occasionally we followed-up targets 
outside our formal bounds, but despite its spectroscopic 
confirmation we cannot include ZwC10628_7_08_0 in the 



rates. 

Second, our control time goes to zero for SN II fainter 
than Mb ~ —14.6 magnitudes, yielding an unphysical 
infinite rate. We exclude these realizations when calcu- 
lating the median rate and its uncertainty, which is effec- 
tively the same as truncating the SNII LF at M ~ —14.6 
magnitudes. This approach is valid because it is inappro- 
priate to calculate the rates of objects to which a survey 
is insensitive. The result is that our rates are for the 
"normal" population of SNe II with Mb < —14.6 magni- 
tudes, and do not include the faint sub-population found 
by Lilla. If we include the realizations of infinite rates 
in the median (e.g. for "All" galaxy types and R200), 
the rate becomes 0.031 instead of 0.026 SNuM with the 
Lilla LF - a difference of < 0.07cr. For the C93a LF, 
the difference is also negligible. In Appendix B we dis- 
cuss a rate calculation method which integrates over the 
SN II LF to avoid the instances of zero control times, and 
explain how it is not appropriate for MENeaCS. 

Third, as discussed in § [2] our sparse light curve sam- 
pling and single epoch spectroscopy means we cannot 
identify SNII subtypes P, L, b, or n. It is likely that the 
subtype distribution in clusters is similar to that in field 
SO-Sbc hosts reported by Lilla, so it is possible that our 
sample contains all SN IIP. For this reason, we also run 
our Monte Carlo using SNII LFs for the plateau subtype 
only. Figure [6] shows that the SNIIP LFs are fainter, 
and that the L, b, and n subtypes populate the magni- 
tude bins brighter than -17.5. This leads to higher rates 
for SNIIP: 0.032 instead of 0.026 SNuM with the Lilla 
LF, which is a difference of < 0.08er. 

4.5. Comparison to Published SN II Rates 

The only previous measure of the cluster SNRu is 
from M08, who compiled five visual and photographic 
z < 0.02 galaxy-targeted SN searches in which "the orig- 
inal galaxy sample was not selected in order to reproduce 
the cosmic average but rather to have a significant num- 
ber of SN detections" . This means that massive galaxies 
were preferentially targeted and associated with galaxy 
clusters later, and that the surveyed cluster mass of M08 
is incomplete. Furthermore, SNe discovered during the 
surveys used by M08 were not all spectroscopically clas- 
sified. MENeaCS, on the other hand, surveyed galaxy 
clusters between 0.05 < z < 0.15 to -R20O7 including the 
stellar mass in faint galaxies and the intracluster popula- 
tion, and uses only spectroscopically confirmed SNe. Al- 
though the MENeaCS survey strategy is better, our sur- 
vey duration is shorter, our discovered number of SNe II 
is smaller, and the Poisson uncertainties are larger. In 
general the M08 rates in Tableland the MENeaCS rates 
in Tabled for the variety of cluster radii and galaxy types 
considered, are consistent at the l-2a level. 

A comparison of the rate of core collapse supernovae 
(SNRcc, which includes Types II and Ibc) between clus- 
ter and field galaxies with similar SFRs can potentially 
reveal an environmental dependence of the initial mass 
function (IMF). For example, consider two galaxies with 
the same total star formation rate: if one has a higher 
SNRcc, then it is forming a larger fraction of > 8 M Q 
stars. We did not discover any SNIbc. The Ibc:II ratio 
is typically 1:2 to 1:4 in cluster and field environments, 
respectively (M08; Li et al. 2011b). Based on this, we 
only expect 1.7-3.5 SNe Ibc, and the probability of ob- 
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serving zero is 3-18%. Thus, we are insensitive to any 
environmental dependence of the IMF at masses > 20 
M Q . A comparison of the MENeaCS SNR H in cluster 
"RS" and "Off-RS" galaxies to the published rates in 
field E/SO and SOa/b galaxies in Table [T] finds a general 
l-2cr agreement. We conclude that SNII cluster rates 
require a more precise measurement in order to firmly 
identify any difference in the cluster IMF. 

5. CLUSTER STAR FORMATION RATES 

We begin our discussion with the first-ever derivation 
of cluster SFR from SNell, which we compare with other 
measurements of cluster SFR, in § 15.11 In § 15.21 we dis- 
cuss the MENeaCS limits on intracluster star formation. 
In § 15.31 we discuss the implications of cluster SF for 
SNela, including whether short-delay SNela contami- 
nate the late-time DTD when measured via the cluster 
SNIa rate; the possibility that all SNela have short de- 
lays; and the potential source of the SNR] a enhancement 
in cluster ellipticals observed by M08. 

5.1. Derivation of Cluster SFR from SN II Rates 

Botticella et al. (2012) present a derivation of the rela- 
tion between the CC SN rate and SFR, which we adapt 
to SNII: 



Fin = K n x SFR. 



(4) 



In the above expression, Ru is the rate of SNII with 
units of SN yr" 1 (whereas SNRjj is in units of SNuM 
= SN (lOV) -1 (lO^M©)" 1 ). The value of R H is cal- 
culated in a similar manner as SNRu, and presented in 
Table El The SFR is in units of M© yr" 1 , and K H is 
the number fraction of all stars formed which explode as 
SNI^SNMq 1 ): 



r u,u (h(rn)dm 



(5) 



The integration limits miji and m u ji correspond to the 
minimum and maximum initial masses of stars which be- 
come SNe II. The lower limit is generally agreed to be ~ 8 
Mq , and evidence is converging towards an upper limit of 
< 20 M (e.g. Smartt ct al. 2009; Dcssart et al. 2010). 
For cf>(m), the initial mass function (IMF), we use the 
Salpeter expression where (j)(rn) oc to 7 , and 7 — —2.35 
(Salpeter 1955). The limits to; and m u are the mass 
range of the IMF, for which we use 0.1 and 100 M©. Un- 
der these assumptions, Ku — 0.0054. By extending to; 
and m u to 0.05 and 200, testing a slightly lower expo- 
nent of 7 = —2.3 (Kroupa 2001), and considering a lower 
muji = 17 Mq, we estimate an uncertainty on Ku of 
±0.001. We account for this with a ~ 18% systematic 
on the cluster SFR. We also calculate the specific star 
formation rate as sSFR = SNRu/Ku, and list them 
in Table El 

We now compare with a selection of previous cluster 
SFR measurements from IR and UV imaging, and opti- 
cal spectroscopy. Chung et al. (2011) used WISE images 
of 72 low redshift (z < 0.1) galaxy clusters at wave- 
lengths 3.4, 4.6, 12, and 22 fim. They found that for 
0.5 < R < 1 -R200, clusters have a mean sSFR ~ 2 
M© yr" 1 (10 12 M©) _1 (see their Figure 3). Our mea- 
surement for "All" galaxy types in clusters, sSFR = 



5.lt3 5 i 8 ±0.9 M© yr" 1 (10 12 M©) _1 , is in agreement with 
this. Yi et al. (2005) used GALEX UV images to look 
for evidence of recent SF in SDSS field early-type galax- 
ies. They found that 1-2% of the stellar mass in ~ 15% 
of bright (Mr < —22 magnitudes) early-types formed 
within the last ~ 1 Gyr. If this is true for clusters, we 
would expect SFR = 2-4 M© yr" 1 in bright red se- 
quence galaxies. If true for fainter red sequence members 
also, the expected total SFR raises to 10-20 M© yr" 1 , 
which is ~ 5 times higher than (but within 2a of) our 
red sequence SFR, 2.1±§ g ± 0.4 M© yr" 1 (lO^M©)" 1 . 
Based on the Ha luminosity from SDSS optical spec- 
tra, Finn et al. (2008) determined that sSFR ~ 9 
M© yr" 1 (10 12 M©)" 1 for star-forming galaxies in low 
redshift clusters. This is in la agreement with our "Off- 
RS" sSFR = 16.6^ 9 6 5 ±3.0 M© yr" 1 (lO^M©)" 1 . 

Finally, we note that our low sSFR in red sequence 
galaxies implies that only ~ 0.01% of the mass in the 
red sequence at z ~ 0.1 is comprised of stars formed 
in the past 50 Myr. This extremely low percentage is 
consistent with the lack of detection in UV and IR of the 
red sequence host of Abcll399_ll_19_0 presented in § [33 

5.2. Intracluster Star Formation 

Direct evidence of intracluster star formation has been 
presented by Sun et al. (2010), who detected a 40 kpc 
long X-ray tail extending from a galaxy in a nearby rich 
cluster. With optical spectra they identified 35 H II re- 
gions along this tail, the furthest of which are 20 kpc 
away from the galaxy (far enough to be defined as in- 
tracluster). Also, simulations of galaxy clusters inves- 
tigating the size and origin of the IC stellar population 
suggest that ~ 30% of the IC stars form a t significant dis- 
tance s from a galaxy dark matter halo ()Puchwein et al.l 
2010). In such simulations most of the IC stars form 
at z > 1, with just a small tail of ~ 1.5% of the final 
IC stellar mass forming during the last ~ 1.8 Gyr (since 
z = 0.15). For the average IC stellar mass of our clusters, 
assuming an IC mass fraction of 16-45% (Gonzales et al. 
2005; Sll), this implies an IC SFR of 5-13 M© yr" 1 , and 
Ru = 0.019-0.049 SNII yr" 1 . This converts into an ex- 
pectation of 0.3-0.8 IC SNe II in the MENeaCS sample, 
which is consistent with our observation of zero. Inter- 
estingly, if the Puchwein et al. (2010) simulations are 
correct, the first detection of an IC SNII would be likely 
in a survey just twice as large as MENeaCS. 

Although we know our upper limit on the IC SFR will 
not be very restrictive, it is the first derived from the non- 
detection of IC SNe II. This is only possible in a complete, 
well characterized survey like MENeaCS. To calculate 
an upper limit for the rate of IC SNe II wit hin R200, 
we us e Poisson statistics for a detection of zero (|Gehrels I 
11986J), and allow fractional values of IC Nu in the Monte 
Carlo rate calculation represented by Equation [T] (with 
Cmc = 1.0). The resulting la upper limit is Ru < 0.15 
SN yr" 1 , which converts to an upper limit of IC SFR < 
28 M© yr" 1 . 

5.3. Implications For Cluster SNela 

In § [1] we described how deducing the late-time SN la 
DTD from the cluster SNRj a as a function of redshift 
must assume that cluster stars all formed in a burst at 
high redshift, and that all cluster SNela have experi- 
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enced a long, > 2 Gyr, delay time. In this section we 
explore the impact of cluster star formation on this as- 
sumption, and its implications for the SNIa DTD. First, 
we discuss our results with respect to the possibility of 
prompt-only SNIa DTD. Next, we estimate the contam- 
ination of "prompt" SNe la to the cluster sample at low 
redshift. Finally, we comment on the implications of our 
cluster SNRu and SFR regarding the source of the en- 
hanced SNRj a in cluster ellipticals presented by M08. 

5.3.1. A DTD of Short Delays Only? 

Maoz et al. (2010) find that the cluster SNR Ia as a 
function of redshift is best and most simply fit by a brief 
burst of star formation at z ~ 3 combined with a SN la 
DTD that peaks at short delays of < 2 Gyr and decreases 
as a power law with a slope of -1 to long delays of ~ 11 
Gyr. This DTD is consistent with theoretical predictions 
of the double degenerate scenario. They also consider a 
DTD of short delays only, congruent with predictions 
for the single-degenerate model. They find it can only 
reproduce the observed cluster SNRj a (z) from < z < 
1.2 if they include ongoing cluster star formation within 
a radius of 1 Mpc. 

Maoz et al. (2010) reject the prompt-only DTD hy- 
pothesis in part because previous surveys found SNe la 
mainly within 1 Mpc and always in elliptical, red se- 
quence galaxies showing no signs of star formation, and 
in part due to other work showing cluster star formation 
was predominantly outside of 1 Mpc. With MENeaCS, 
we have shown that SNe la do occur outside of 1 Mpc, 
and in blue cluster galaxies (Figure [3J and S12). We have 
also shown that SNe II, and therefore star formation, oc- 
curs inside of 1 Mpc and in red sequence galaxies. 

Qualitatively, this suggests it would be premature to 
rule out a DTD of short delays only, and the single- 
degenerate model as the sole scenario, based on observed 
SNe la cluster rates. Quantitatively, Maoz et al. (2010) 
find that a SNIa DTD of short delays only requires a 
rate of ongoing star formation SFR ~ 175 M Q yr -1 in 
the central regions of galaxy clusters. This is ~ 3 times 
higher than our cluster SFR inferred from SNe II, but our 
uncertainties on cluster SFR are large, and 175 M Q yr -1 
is actually just within the la upper limit. Ultimately, 
although our cluster SFR rate is low, we cannot rule out 
a DTD of short delays only. 

5.3.2. Fraction of "Prompt" Cluster SNe la 

The SNIa rate per unit mass, SNRi a , is a convolution 
of star formation history of the surveyed galaxy sample 
and the SNIa DTD. As an oversimplified parametriza- 
tion, it can be expressed as the sum of "delayed" and 
"prompt" components represented by c onstants A and 
-B, where SNR Tn = A + B x sSFR (iMannucci et al.1 
[20051: iScannapieco fc Bildste n 2005; Sullivan et al.ll2006l) . 
Based on this, the fraction of SNe la expected to be 
associated with the prompt component is: /prompt = 
(B x sSFR)/SNR Ia . We use the B value from Sullivan 
et al. (2006), 3.9 ± 0.7 x 10~ 4 SN yr" 1 (Moyr" 1 )" 1 . 
For MENeaCS red sequence cluster galaxies, our de- 
rived sSFR is in Table E| and SNR Ia = 0.04lt°;^+°-0°5 
SNuM (S12). This reveals f prompt = 0.02±g;g| in cluster 

red sequence galaxies; similarly, f pr0 mpt = 0.05^ O O5 for 
"All" cluster galaxies. 



While this minimal ^2% contamination indicates that 
all of the MENeaCS SNe la likely experienced a long de- 
lay time, we cannot rule out that up to 7% of SNe la in 
red sequence hosts exploded with short delay times (the 
la confidence level). If so, then the "RS" SNRi a for the 
delayed SNe la only would be ~ 0.038 instead of 0.041 
SNuM. This is a difference of 0.003 SNuM or, expressed 
in terms of the uncertainty on SNRi a , ~ 0.2ct. Ulti- 
mately, we find the potential maximum contamination 
from "prompt" SNe la is less than our statistical uncer- 
tainty on the cluster SNIa rate. However, this may not 
be the case for higher redshift SNIa cluster surveys. 

5.3.3. The Enhanced SNRi a in Cluster Ellipticals 

The rate of SNe la in cluster early-type galaxies was 
found to be a factor of three higher than the rate in field 
early-types by M08: SNR Ia = 0.066 SNuM compared to 
0.019 SNuM, respectively. They report that if this excess 
is from the "prompt" component, they would only expect 
~ 2 SNe II - consistent with their detection of no SNe II 
in E/S0 hosts. In S12, we report the cluster red sequence 
SNRia = 0.04lt° °^+°-°°5 SNuM, which is an excess of 
0.022 SNuM over the rate in field early-types from M08 
(but given the uncertainties still consistent). Assuming 
this is entirely from the prompt component, and using 
the B value from Sullivan et al. (2006), we find an im- 
plied sSFR = 56±% M Q yr" 1 (IO^Mq)" 1 in cluster red 
sequence galaxies. This is a factor of ~ 28 times greater 
than our observed SFR in Table (5[ If we instead con- 
sider our "All" galaxy sample the implied sSFR = 56^34 
M Q yr -1 (10 12 M Q ) _1 is a factor of ~ 11 times greater 
than the observed SFR in clusters. Given that the rela- 
tion between SFR and SNRu is direct (Equation [4}, 
then we should have observed an order of magnitude 
more SNe II in MENeaCS. Our data suggests that either 
the SN la rate enhancement in cluster ellipticals does not 
exist, or it is not due to recent star formation. 

6. CONCLUSION 

In this paper we present the 7 SNe II discovered in 
0.05 < z < 0.15 rich galaxy clusters by MENeaCS. Our 
sample also includes one SNII in a red sequence galaxy 
which shows no clear evidence of recent star formation 
in its multi-wavelength properties. This illustrates the 
danger of using host morphology to classify SNe in lieu of 
expensive spectroscopy time. The simplest explanation 
is that undetectable levels of star formation exist in the 
elliptical host. If that is not the case, it leaves open the 
possibility of a rare other channel to SNe II with a long 
delay time. 

With the MENeaCS sample we make the first mea- 
surement of SNRu from a survey which is both cluster- 
targeted and complete to i?200- We also make the first 
derivation of cluster SFR from SNRu, and find that it 
agrees with SFR measurements for cluster galaxies and 
field ellipticals from IR and UV photometry, and Ha line 
emission. We show how these low levels of cluster star 
formation imply that a small fraction of cluster SNe la 
may have experienced a short delay time. However, we 
find their influence on the cluster SNRi a is within sta- 
tistical uncertainties, and does not undermine the use of 
low redshift cluster SNRi a (z) to derive the delay time 
distribution for SNe la. 
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TABLE 1 
Published SNII Rates 



Environment 


Galaxy 
Types 


Rate 
SNuB a 


Rate 
SNuM b 


Mannucci et al. (2008) 


Cluster 

Cluster 

Cluster 

Cluster 

Field 

Field 

Field 

Field 


Total 

E/SO 

SOa/b 

Sbc/d 

Total 

E/SO 

SOa/b 

Sbc/d 


n 9 o+U.09 
u - zo -0.07 

< 0.071 

o-87l2;S 

44+ - 08 

< 0.08 

°- 36 -o.io 
0.83t°;£ 


U.U(Z_g Q21 

< 0.017 

061+ 068 

U.UDl_g 03 g 

610+ ' 290 

U.DlU_g jog 

174+0.032 
o-J-I"±_q 027 

< 0.020 

130+ 049 

u.±ou_ 037 

652+° 164 

u - od -0.134 


Li ct al. (2011b) 


Field 
Field 
Field 
Field 
Field 
Field 
Field 


E 

SO 

Sab 

Sb 

Sbc 

Sc 

Scd 


< 0.014 
0.200±o °g 5 (0.006) 
0.266l2;2tl (0.098) 
0.282t C o.°)37 (0.106) 
0.466 + 2'EJs2 (0.134) 

nfi4 q+0.088 ,+0.364-. 

n 7QK+0.097 ,+0.386-, 
u -' ao -0.086 V-0.135-* 


< 0.003 
0.005±°;°°* (0.001) 
0.098l2;21| (0.035) 
0.144±^°lg (0.055) 
0.355l°;g|l (0.098) 

n K47+0.075 ,+0.245-, 

+0.116 ,+0.342-, 
U-'D'_0.102 ^.-0.154/ 



a SNuB = SNe(100 yr 10 10 is,©) -1 
b SNuM = SNc(100 yr 10 10 M©)" 1 



TABLE 2 

MENeaCS Type II Cluster Supernovae Spectroscopy Summary 



MENEACS ID 



UT Date 



Telescope/ 
Instrument 



Cluster Galaxy 



SNID 



SNID Template 
Type, Name 



SNID 
Phase (<t) 
days 



Exposure 

Time 
seconds 



Abellll9_5_24_0 

Abell 1795.8.08.1 

Abell399_ll_19.0 

Abclll651.7.05.3 

ZwC10628.7.08.0 

Abcll2443_5.19.0 

Abell990_6.13.0 



2009-09-19.29 
2009-06-15.27 
2009-12-19.27 
2009-06-15.20 
2009-09-19.49 
2009-06-15.42 
2009-03-16.11 



MMT/BCS 
MMT/BCS 
MMT/Hecto 
MMT/BCS 
MMT/BCS 
MMT/BCS 
MMT/Hecto 



0.044 
0.063 
0.072 
0.085 
0.081 
0.108 
0.144 



0.0480 
0.0626 
0.072 



0.1106 
0.1425 



0.050 (0.005) 
0.0636 (0.0054) 
0.0676 (0.0025) 
0.0739 (0.0050) 
0.0759 (0.0050) 
0.1053 (0.0056) 
0.1422 (0.0039) 



IIP, SN92H 
IIP, SN04ct 
IIP, SN04ct 
IIP, SN04et 
IIP, SN04et 
IIP, SN99em 
IIP, SN04et 



39.5 (88.0) 
21.4 (55.8) 

60.1 (19.2) 
26.9 (20.9) 

48.2 (37.4) 
24.0 (15.5) 

14.3 (10.8) 



1200.0 
900.0 
3600.0 
900.0 
900.0 
2400.0 
3600.0 



TABLE 3 

MENeaCS Type II Cluster Supernovae Host Properties 



MENEACS ID 


a 


S 


r 


g-r 


A(g - r) RS 


-R-clus 


-f^clus 




J2000.0 


J2000.0 


mag 


mag 


mag 


kpc 


R200 


Abcllll9_5.24_0 


00:55:39.69 


-00:52:35.9 


16.53 + 0.02 


0.40 ± 0.04 


-0.41 


1280 


0.78 


Abelll795.8.08.1 


13:48:38.56 


+26:22:18.5 


20.23 ± 0.02 


0.31 + 0.04 


-0.37 


980 


0.46 


Abell399.11.19.0 


02:57:16.58 


+13:08:34.3 


16.28 + 0.02 


0.85 ± 0.04 


-0.02 


920 


0.49 


Abclll651.7.05.3 


12:57:23.39 


-04:33:52.7 


22.26 + 0.06 


0.22 + 0.10 


-0.45 


3530 


1.73 


ZwC10628.7.08.0 


06:31:03.31 


+24:49:17.8 


16.11 + 0.02 


0.79 ± 0.05 


-0.15 


1170 


0.73 


Abcll2443_5.19.0 


22:25:07.59 


+17:35:00.9 


17.35 ± 0.02 


0.56 ± 0.04 


-0.43 


2350 


1.47 


Abell990.6.13.0 


10:23:18.58 


+49:05:17.8 


20.97 + 0.03 


0.32 ± 0.05 


-0.63 


730 


0.39 
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TABLE 4 
MENeaCS Cluster SN II Rates 



Galaxy Radius Ngjvj SNII Rates 

Types SNuB a SNuM b 

With the SNII LF of Li et al. (2010). 

"J:^^) 0.032^;^(stat)l ;^(sys) 

3±8lg|? (stat)l g : «3 (sys) 0.0261°;°!! (stat)l g;°S° (sys) 

^°;^ (stat) +o.oo5 (sys) o.on+o.035 (stat) +o.ooi (sys) 

n^( S tat)l°;°° 3 (sys) 0.007i°;°J 7 (^at)l°;°°?(sys) 

iS:^(stat)l°;°«(sys) 0.129i°;t«(stat)l°;°f 7 (sys) 

+° jf (stat)l °;gf 9 (sys) 0.0801 g;^ (stat)±°;8H (sys) 



All 


1 Mpc 


3 


0.137 


All 


R200 


4 


0.093; 


RS 


1 Mpc 


1 


0.044 


RS 


R200 


1 


0.028 


Off RS 


1 Mpc 


2 


0.346 


Off RS 


R200 


3 


0.220 



With the SNII LF of Cappcllaro ct al. (1993aJ~ 



All 


1 Mpc 


3 


0.112 


All 


R200 


4 


0.088 


RS 


1 Mpc 


1 


0.035; 


RS 


R200 


1 


0.023 


Off RS 


1 Mpc 


2 


0.305 


Off RS 


R200 


3 


0.194 



lo 073 («tat)l!-- (sys) 0.035 ;. ;-- (stat)l^ (sys) 

1 063 ( stat )± 0T1 ( s y g ) 0.026l°-«8 (stat)l°'°J° (sys) 

+0.165^, ,\+0.004/ \ nnno+0.041/. ,.,+0.001/ \ 

-0.035 ( stat )-o. 003 (sys) 0.008_ 0()8 (stat)_ 001 (sys) 

±C:«23 (rtat)l g;g*» (sys) 0.005±°;g£ (stat)lg;gJ 2 Q (sys) 

tlfol (stat)l°;°^ (sys) 0. 106ig;?& (stat)l g;gj» (sys) 

±8?|1 (stat)lg- 2 ^ (sys) 0.0671°/^ (stat)±g-g£ (sys) 



SNuB = SNc(100 yr 10 10 L s , )~ 
' SNuM = SNe(100 yr 10 10 M )" 



TABLE 5 

MENeaCS Cluster Star Formation Rates 



Galaxy 
Types 


N S N 


Rii 
SN yr" 1 


SFR 

M yr- 1 


sSFR 

M yr" 1 (IO^Mq)- 1 


All 

RS 

Off RS 


4 
1 

3 


0.308l° ;?- 
0.072lg;Jf 2 
0.225lg;f 3 6 


56.6lj3 7 7 3 (stat) ± 10.2(sys) 
13.3lg 8 g 4 (stat) ± 2.4(sys) 
4L3 -24.i( stat ) ±7.4(sys) 


5.ll i 3 5 1 8 (stat) ±0.9(sys) 
2.0lo g(stat) ±0.4(sys) 
16.6lg 9 6 5 (stat) ±3.0(sys) 
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Fig. 1. — Spectra for our cluster SNell in black, with the best fitting SNII template spectrum from SNID in red (see Tabled- 
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Fig. 2. — The offset of each SNII (black) and SNIa (gray) host galaxy from its cluster's red sequence, A(g — r)ns, versus the host galaxy 
r-band apparent magnitude (left) and the clustercentric offset in units of i?200 (right). Dashed lines represent the median scatter in the 
red sequence over the whole MENeaCS sample. 




Host Color (g— r) Clustercentric Distance (R/R Z oo) 

Fig. 3. — Left: The normalized cumulative distributions of SNIa (gray) and SNII (black) host g — r colors, for SNe within -R200- A two- 
sided KS test finds there a < 10% probability that these are drawn from the same underlying distribution. Right: normalized, cumulative 
g-band luminosity in red sequence galaxies (dashed), and the cumulative number fractions of SNe la (gray) and SNe II (black). Within 
^200: the distribution of projected clustercentric offsets for SNe la appears to "follow the light" in the red sequence better than SNe II, as 
expected. 




Fig. 4. — Left: The host of SNII Abell399_ll_19_0, with the SN marked as a red cross, and smoothed isophotal contours to highlight this 
galaxy's elliptical morphology. Right: A partial spectrum of the red sequence host galaxy of SNII Abell399_ll_19_0, showing the relevant 
emission lines. Line intensities are estimated by integrating the flux over the line width. 
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Fig. 5. — The best fitting template from Superfit (red) to the observed spectrum of Abell399_ll_19_0 (blue, with best fittin g host spectra 
subtracted), when considering only SNII (left) and only SNIa (right) with the spectral template matching code Superfit l|Howell et al.l 
[20051) . The best fit SNII template is a much better match than the best SN la template, in particular to the P-Cygni profile at ~ 7000 A. 
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Fig. 6. — SNII luminosity functions. Black, the volume-limited SNII luminosity function observed by Li et al. (2010), for host galaxies 
of type SO-Sbc. Blue, the SNII luminosity function used by Cappellaro et al. (1993a), normalized to the same total number as the Li et 
al. (2010). 
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A. DERIVATION OF THE INCLINATION CORRECTION 

A lower SNRjj has been observed in high inclination spiral galaxies relative to low inclination spirals (e.g. Cappellaro 
et al. 1993b; Cappellaro et al. 1999). This difference is attributed to dust obscuring SNe when our line of site goes 
through the disk. It was suggested by Cappellaro et al. (1993b) that this only affects photographic surveys. However, 
Li et al. (2011b) is a CCD survey and they report a strong trend between SNRu and inclination for late-type spirals, 
but a very mild trend for early-type spirals. 

Although the SNII LF of Lilla was not corrected for host extinction, using it does not mean we have included an 
inclination correction factor. For galaxy targeted SN surveys the inclination of each spiral galaxy is inferred from its 
shape, and the inclination correction derived from the observed rates. However, without redshifts for all galaxies in 
the CFHT MegaCam field of view, we cannot identify which spirals belong to the cluster and our inclination correction 
must be statistically estimated. 

To derive our inclination correction factor, we first describe the simple case of having two samples of galaxies, 
"regular" and "inclined". The fraction of SNII host galaxies which are "regular" is fn(r), the fraction which are 
"inclined" is an d the fraction of all occurring SNe II that are detected is fsN d (r) and fsN d (i)- The inclination- 

corrected number of SNe II, Nn tC , given the total number observed, iV/j, is Nu tC = Ci nc Nn where the inclination 
correction factor Ci nc is: 

Ci nc = fH(r)-? 1" fH(i)-? • (6) 

JSN d (r) ' JSN d (i) 

Instead of simply "regular" and "inclined" , we consider 7 different galaxy types: ellipticals, and spirals of early (Sab) 
or late (Scd) types which are face-on (f, < i < 40), inclined (i, 40 < i < 75), or edge-on (e, 75 < i < 90). From 
our observations we know the approximate fractions of SNe II hosted by elliptical and spiral galaxies are /he ~ 1/7 
and fns ~ 6/7. Based on the ratio between LOSS rates in field Sab and Scd spiral galaxies listed in Table [TJ SNe II 
occur ~ 10 times more often in Scd than Sab spirals. We break fns down into components fnsab — 0.086 and 
fused = 0.771. Maltby et al. (2011) find he number fraction of cluster spirals which are face-on, inclined, and 
edge-on is 0.16, 0.66, and 0.18 respectively. The fractions of SNII hosts in each of our 7 considered galaxy types 
are: f HE = 0.143, fusabtf) = 0.014, f H Sab(i) = 0.057, f H Sab(e) = 0.015, f H Scd{f) = 0.0123, f H Scd{^ = 0.509, and 

fHScd(r) = 0.139. 

Assuming all SNe II in face-on spirals are detected, we derive the fraction of SNe II detected in inclined and edge- 
on spirals from the rates presented in Table 3 of Li et al. (2011b). MENeaCS and LOSS are both multi-epoch 
CCD SN surveys at low redshift, and it is reasonable to expect a similar inclination effect. LOSS did not target 
galaxies expected to produce extremely extincted transients such as very dusty starburst galaxies, but neither does 
MENeaCS since starbursts are not generally found in rich galaxy clusters. They consider the same 7 galaxy types 
as described above, and from their work we find fsN d (E)=l-0, /sjv d (s o 6(/)l=l-0, fsN d {Sab(i))=0 -78, fsN d (Sab(e))=0-6, 
/sjv <j (S<j6(/))=1-0, fsN d (Sab(i))=0-63, /sAf <J (Sab( e ))=0.32. This results in a MENeaCS inclination correction factor of 
Cind— 1.62. By varying the component terms between their minimum and maximum values (e.g. the LOSS rates 
uncertainties), we estimate a ~ 5% uncertainty on Ci nc and add this to our systematic uncertainties. Note that Ci nc 
does not apply to our rates in red sequence galaxies. 

B. APPROXIMATION TO A RATE INCORPORATING AN "EFFECTIVE" CONTROL TIME 

The occurrence of unphysical rates can be avoided if, instead of choosing randomly from the SNII LF for each 
realization of the Monte Carlo, an effective control time is calculated using a sum over the components of the LF. In 
this method, Equation [1] becomes: 

where /, is the weight of the i th bin of the SNII LF, Nzii is the number of observed SNe II in the i th bin, and Atij 
is the control time for SNII from the i th bin. This method is derived in the Appendix of Leaman et al. (2011), and 
is used by LOSS. Since MENeaCS light curves are sparsely sampled, we do not have Nuj, and this method is not 
appropriate for our dataset. We implement an approximation to this method by taking iVjj outside the sum over i, 
and bringing Y^ifi inside . to calculate an "effective" control time, integrated over the SNII LF, for every epoch. 
This style of rates calculation is more processor intensive, and for this trial we only run 100 realizations. The resulting 
rate for "All" galaxy types within i? 2 oo is SNRu = 0.023^0 ^5 SNuM, and no unphysical rates are encountered during 
the Monte Carlo. The upper uncertainties are smaller than those for our main method, SNRu = 0.026±Hi8> 

because 

the SN II LF does not contribute to the width of the rates distribution from the Monte Carlo realizations. We reiterate 
the subtle point that this is not appropriate for a survey such as MENeaCS in which the observed LF, Nn,i, is not 
well characterized. 



